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We study a quantum quench in a system of two coupled one-dimensional tubes of interacting 
atoms. After the quench the system is out of equilibrium and oscillates between the tubes with a 
frequency determined by microscopic parameters. Despite the high energy at which the system is 
prepared we find an emergent long time scale responsible for the dephasing of the oscillations and a 
transition at which this time scale diverges. We show that the universal properties of the dephasing 
and the transition arise from an infrared orthogonality catastrophe. Furthermore, we show how this 
universal behavior is realized in a realistic model of fermions with attractive interactions. 
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Ultracold atomic systems provide a unique laboratory 
for investigating noncquilibrium dynamics in correlated 
quantum matter. A common experimental approach is 
to prepare the system in a simple known quantum state, 
which is not an eigenstate of the system Hamiltonian, 
and observe the ensuing time evolution [H, 0, El- This 
is known as a quantum quench. The near perfect isola- 
tion of the atoms from the environment allows to con- 
centrate on the intrinsic many-body dynamics following 
the quench. The focus of this Letter is on the question 
whether such dynamics can give rise to emergent struc- 
tures or universal dynamical modes. 

The theoretical description of quench phenomena poses 
a fundamental challenge. Generically, the system is in- 
jected with a high energy density. Consequently, the time 
evolution involves all energy scales. A natural expecta- 
tion in this case is that the subsequent dynamics would 
be highly complex and nonuniversal: completely domi- 
nated by the microscopic scales at short times; at longer 
times too complicated to capture in a reasonable theory. 

Yet there is accumulating evidence to the contrary. 
First, emergent phenomena may certainly arise when the 
energy of the initial state is not too large. In such cases a 
low energy effective theory appears to give a correct de- 
scription of the quench and bring out universal aspects 
of the dynamics. A good example is the decoherence 
dynamics of a pair of one-dimensional condensates,pre- 
pared in a state with well defined relative phase [2j, |j| . 
The dynamics in this case can be understood rather sim- 
ply within the universal Luttinger liquid description of 
the condensates 0, @ . The same analysis holds for the 
time evolution of coherence between two internal spin 
states in a one dimensional condensate starting from a 
perfectly coherent state Q. Related theoretical ques- 
tions of quench dynamics within conformal field theories, 
of which the Luttinger liquid is a particular example, are 
discussed in Refs. @, 

Second, and more surprisingly, there is numerical evi- 
dence for the existence of emergent dynamical modes of 
very long and possibly diverging time scales, even when 



the initial state in the quench is at high energy. Such 
behavior was seen, for example, in the time evolution 
of the staggered mag netization in a quantum quench of 
a spin- 1/2 chain [lfj. In a wide regime the staggered 
moment displays rapid oscillations, which clearly reflect 
the microscopic magnetic exchange interactions. On the 
other hand, the decay of the oscillations occurs on a much 
longer time scale, which even diverges on approaching the 
XY limit of the spin chain. Such behavior is highly sug- 
gestive of a theoretical description, which separates out 
the fast dynamical modes in order to focus on the emer- 
gent long time dynamics. We do not know, at present, 
of a consistent way to accomplish such separation in the 
spin model of Ref . [n| ■ 

In this Letter we achieve separation of scales within a 
related model of a quantum quench, which is relevant to 
experiments with one dimensional bosons or fermions in 
coupled tubes, Fig. [1] Consider a pair of tubes, where 
in the initial state one of the tubes is at a much lower 
potential. The ground state corresponds to the quan- 
tum liquid filling that tube. How does this ground state 
evolve when the empty tube is suddenly brought to near 
resonance with the filled tube? 

The regime of interest to us in this work is when the 
tunnel coupling between the two tubes J, defines the 
largest frequency scale in the problem. That is, J is 
much larger than the chemical potential fi. In this regime 




FIG. 1: Two tubes detuned by A and coupled via a tunneling 
amplitude J. For |A/J| ^> 1 the ground state is given by 
all particles residing in the lower tube. We investigate the 
behavior of the density after a sudden ramp to A « 0. 
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the problem defines a direct many-body generalization of 
the usual Rabi oscillations of a single particle in a double 
well. The microscopic Hamiltonian of the system is given 
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Here, t q are Pauli matrices acting in the pseudo-spin 
"tube-space" and n = L,R (left/right) is the tube index. 
A is the detuning between the tubes after the quench. 
The Hamiltonian is written for fermions with mass m and 
spin it =t, 4; Pna is the local density operator. For bosons 
the Hamiltonian takes an identical form only dropping 
the spin subscript. The effective contact interaction g 
depends on the three-dimensional scattering length a s 



and the transverse confinement 

In absence of interactions, under the influence of Ho, 
all particles perform Rabi oscillations of frequency £1 — 
h^ 1 ^/ J 2 + A 2 /4, perpetually in phase with each other. 
In general, when such rotating spins are coupled to a 
dynamic environment, the fast time-scale ft^ 1 is supple- 
mented by slow time-scales Tu 2 ), responsible for relax- 
ation and dephasing of the Rabi oscillations, respectively. 
Here, these slow time scales emerge due to an intrin- 
sic mechanism. The fluctuations of the one dimensional 
quantum liquid will supply the sought slow dynamics in 
a tractable way. 

Before proceeding we outline the general strategy for 
solving this problem and our main results. In order to 
separate out the rapid oscillations, we work within a ref- 
erence frame rotating with the pseudo-spins (interaction 
picture with respect to Hq). In this frame the two modes 
corresponding to left and right tubes are replaced by the 
filled, in-phase rotating mode (mode 1) and the orthog- 
onal, out-of-phase rotating mode (mode 2). Without in- 
teractions all particles rotate perpetually in mode 1. The 
interactions give rise both to an effective detuning A ctf , 
as well as particle hopping terms between the two modes. 
The particle current from the in-phase rotating mode to 
the out-of-phase mode is the dephasing rate. 

For small effective detuning we can compute the de- 
phasing rate by treating the quantum liquid in the filled 
tube, within its low energy effective theory, the Luttinger 
liquid. The dephasing rate as a function of effective de- 
tuning is found to be a power law T% ~ ( A off / p) 13 , where 
the exponent /3 can be positive or negative and reflects 
the critical fluctuations of the Luttinger liquid. A transi- 
tion to positive (3 corresponds to an orthogonality catas- 
trophe [13], whereby tunneling of particles out of the Lut- 
tinger liquid (dephasing) is suppressed at zero detuning 



by the critical fluctuations. For negative (3 the tunnel- 
ing ( dephasing) rate is enhanced rather than suppressed 
by the fluctuations. We show that in an experimentally 
relevant system of fermions with attractive interactions, 
the exponent (3 can be tuned from negative to positive 
by tuning, for example, a s . 

We now turn to derive the main results. The transfor- 
mation to the rotating frame is accomplished by 
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In this frame Ho is eliminated and the Hamiltonian is 



H = H 1 +H 2 + H 12 {t). 



(5) 



The intramode Hamiltonians H\ and H 2 are of the same 
form as in Eq. |T]). The coupling between the 

two modes Hi 2 (t) arises because of the interactions and 
contains three different contributions: a single particle 
hopping h lp ~ ga{t)J2a ^L^aiPi-v ~ P2-*) + H - c -; 
a pair hopping term h 2p ~ gb(t)ipl^tp\^ip 2 ^tp 2 ^ + H.c; 
and an inter-mode interaction h int ~ gc(t) p\ a p 2 -o + 
Wio^i-o^X-o^ho- Note that if the interactions were 
symmetric under SU(2) rotations in the tube space then 
H\ 2 would vanish. However, the fact that particles in- 
teract only within the same tube constitutes a strong 
breaking of this SU(2) symmetry. 

The coupling constants appearing in H\ 2 (t) are time 
dependent with a periodicity corresponding to a fraction 
of the Rabi period T. Since the Rabi frequency is 
the highest frequency scale in the system we can treat 
this time dependence within a systematic expansion in 
p/fl [14l ? ]. This approach is the quantum equivalent of 
the Kapitza pendulum problem [l5|. The leading order 
contribution is obtained by replacing a(t), b(t), and c(t) 
by their average value 




where a = A/ J. The oscillatory time dependence of the 
interaction leads to sub-leading corrections at order /i/Sl 
to the calculation of T 2 . 

We focus on the effective time independent problem 
relevant for calculation of T 2 . It consists of a full tube cor- 
responding to the particles in the in-phase mode |1) and 
an empty tube representing the complementary mode |2). 
For the purpose of computing the tunneling rate, the 
empty mode can be represented by 
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Here pT 2 = c int gpo/2 is a Hartree shift due to the effective 
interaction between particles in the two modes. 
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The full mode with density po experiences an interac- 
tion induced shift in its chemical potential p,"" as well. 
We discuss the exact form of this shift away from the 
Fermi energy ep = it 2 h 2 p\ / '8m below. Within a micro- 
scopic model we will also show that tuning the interac- 
tions can lead to an effective resonance between the two 
modes. Namely A off = p^ ~ Mi" = 0- 

For the formulation of the long-wavelength Luttinger 
liquid theory near the effective resonance we intro- 
duce the slow fermion modes ip±, such that ipi,j(x) ~ 
ijj +t(T {x)e lkFX + ip_^(x)e~ lkFX , and kp is the Fermi mo- 
mentum. As usual we express the Hamiltonian ([5]) using 
these fields and neglect terms which oscillate as e mkpX . 
In particular, for this reason the single particle term 
oc a int disappears from the long-wavelength theory. 

Next, we bosonize the slow fields via ip±,<j{ x ) K 
exp{-i[±9 p - ip p + a(±9 a - <Pa)]/V2} M- We shall 
concentrate on the case of attractive interactions, where 
fermions form either (singlet) Cooper pairs or tightly 
bound molecules [13] • From the point of view of the long- 
wavelength theory, the fermions form a Luther-Emery 
liquid with a spin gap A^ [l8l ]. The low energy degrees 
of freedom therefore lie in the charge sector. 



H? = ^rj dx \l< p [d xVp {x)] 2 + ^[d x e p (x)f 



(7) 



The phonon velocity v p is related to the Fermi velocity 
vp = Trhp Q /2m by K p — vp/v p ; later on we derive an 
explicit expression for K p from the microscopic theory. 
The bosonized form of the tunneling terms between the 
two modes is given by 



H\2 = —2g / dx I c, 



i(V26 



cos 
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The set of Eqs. ©-((H]) allows us to calculate the evolu- 
tion of the Rabi oscillations after the quench. We remark 
that due to the spin gap A CT , the first term in ([5]) is ir- 
relevant for A cff < A(j . In the second term of (J8J) , the 
spin-part of the coupling leads to a logarithmic correc- 
tion of the "coupling constant" gb int [16J. Since we are 
only interested in the leading behavior of the dephasing 
rate, we can focus purely on the charge sector. To second 
order in the coupling gb iat , we obtain via a Fermi golden 
rule calculation (FGR) the dephasing rate 

1/T% ~ / dkdqdujdi> A2(k,uj)A2(q,i l ) 

x Ai(—k — q, —u) — v). (9) 
The spectral densities entering the FGR expression are 
A 2 (k,uj) = 2ttS(lj - h 2 k 2 /2m- p), (10) 



Ax{k,oj) = — 



1 



ep \uj 2 — h 2 v 2 k 2 



(11) 




FIG. 2: "Phase diagram" of dephasing rates, (a) Dephasing 
rate for the coherent Rabi oscillations after a sudden quench 
to a fixed effective detuning A off (arb. units). For 7 < 
(BCS) the dephasing rate is suppressed on the effective reso- 
nance, while for 7 > (BEC) it diverges, (b) Phase diagram 
in the plane of I/7 — A/J. To the left is the BEC regime, 
where the dephasing rate diverges, and to the right the BCS 
regime where it is suppressed. The locus of effective resonance 
(A ctf = 0) is marked by the lines enclosing the white area, (c) 
The Luttinger parameter as a function of I/7 [17]]. 



Performing the integral in ([9]) we obtain the dephasing 
rate, which is one of our main results 



l/T 2 ~%^(A eIf /e F )- 1+1 



(12) 



The most striking aspect of this result is that the point 
K p = 2 represents a transition in the nonequilibrium dy- 
namics. While for K p > 2, the Rabi-oscillations are over- 
damped, for K p < 2 the dephasing rate is suppressed as a 
power-law as A cff — > 0. That is, the time scale X2 diverges 
[see Fig. Ufa)]. The universal behavior of T2 is generated 
by an infra-red singularity, or orthogonality catastrophe, 
due to coupling with critical charge fluctuations in the 
Luttinger liquid. The effect is analogous to blocking of 
tunneling through a point contact between one dimen- 
sional electron liquids Scaling analysis, similar to 
Ref. [19(| , shows that coupling between the two rotating 
modes is irrelevant at all orders. Therefore absence of 
dephasing in the limit A off — » is valid beyond the FGR 
analysis. 

All this is strictly true only within our static (time av- 
erage) approximation. The neglected rapid oscillation of 
the coupling in the rotating frame will inevitably lead to 
some dephasing, as will inhomogeneities and finite tem- 
perature. In the opposite regime of K > 2, the divergence 
of the dephasing rate signals failure of FGR as A cff — > 
and it will be cut off in reality. The practical implications 
however remain unchanged: overdamped oscillations for 
K > 2 versus strongly suppressed dephasing for K < 2 
as A off — > 0. 

It is interesting to contrast the dephasing of oscilla- 
tions in the interacting quantum liquid with the usual 
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problem of decoupled two-level systems connected to a 
reservoir. In both cases "T2" processes correspond to 
flipping a two-level system from in-phase to out-of-phasc 
precession. Contrary to the usual case, here the in- and 
out-of-phase modes arc not strictly degenerate. Trans- 
ferring a particle into the out of phase mode can involve 
a change in the interaction or axial kinetic energy. The 
excess energy is absorbed by excitations of the quantum 
liquid. 

In the following we show how a low energy theory with 
K p G [1, 00], A s > 0, and A eft — > is realized for fermions 
confined to quasi-one-dimension. We follow the analysis 
of Fuchs et al. [17j . who discussed an exactly solvable 
model relevant to the BCS to BEC cross-over of quasi 
one-dimensional fermions. The model is parametrized 
by a dimcnsionlcss interaction constant 7 = g/(fi 2 po), 
where g is the coupling constant of the Hamiltonian fTJ) . 
For I/7 — > —00 the system is in the weakly interact- 
ing BCS regime with a spin gap A a oc exp[— 7r 2 /2|'y|] . 
For I/7 > the Hamiltonian (TTJ) docs not support a 
bound state. However, this is an artifact of the strictly 
single channel description. In reality the bound state, 
and the spin gap go up to the transverse trap frequency 
as I/7 — > + violating the single channel assumption, 
this is the so called confinement induced resonance (CIR) 
[TlL Hit ]. Fuchs et al. therefore considered the model 
(TT]) supplemented with a bound state at ej, = — 00 (or 
A CT = 00) for 7 > 0. We can extract the relevant pa- 
rameters for the low energy theory from their exact so- 
lution. In the crossover regime 1/7 « 0, the Luttinger 
parameter K~ s» (1 — l/j)/2 and the chemical poten- 
tial (if « ep(l/4 — I/37). It is interesting to note that 
the transition between over-damped Rabi oscillations to 
suppressed dephasing for K p < 2 occurs exactly at the 
CIR point 1/7 = 0, see Fig.^c). 

To find the location of the effective resonance A off = 0, 
we compute the interaction-induced detuning 



A off = - (if « e F 



1 

3H 



7Cint — 



(13) 



We see that A off = requires 1/7 < 1, which is in the 
strongly interacting, or cross-over regime. The need for 
strong interactions is not surprising because the inter- 
action effects have to overcome the Fermi energy of the 
full mode. In Fig. [DJb) we show the locus of the effec- 
tive resonance on the two sides of the "transition" in 
the plane of bare detuning versus interaction parameter. 
The very different trends on the two sides of the BCS to 
BEC crossover can be observed even if one does not tune 
exactly to A cff = 0. The white region in Fig. [D(b) corre- 
sponds to A cff > where within our FGR calculation no 
particles are transferred out of the full mode. 

We now turn to briefly discuss a similar setup of spin- 
less bosons with repulsive interactions. As mentioned, 
this situation is described by Eq. fT]) dropping the spin 
index. The transformation to a rotating frame remains 



unchanged. Because of the absence of a Fermi energy 
the effective detuning is different, and in the Hartree ap- 
proximation given by A 8 « = pog[(l — |c jDt |) — 2|c int |]. The 
low energy theory of the interacting bosons is also a Lut- 
tinger liquid. However, in contrast to the fermion case, 
the single particle tunneling does not vanish on going to 
the long wavelength limit. In fact it is the most relevant 

term leading to I/T2 ~ A rff 2K . Here, K is the Lut- 
tinger parameter controlling the decay of correlations in 
the filled mode. An orthogonality catastrophe, similar to 
the Fermions, occurs for K < 1/2, which requires bosons 
with longer range (e.g. dipolar) interactions. Alterna- 
tively in the presence of a commensurate lattice poten- 
tial the transition to a gapped Mott insulator would be 
accompanied by a transition to suppressed dephasing. 

In conclusion, ultracold atoms oscillating between a 
pair of one dimensional traps following a quantum quench 
constitute a natural many-body generalization of Rabi 
oscillations. We described universal behavior in the de- 
phasing rate seen as power-law dependence on an effec- 
tive detuning parameter. A transition from enhanced to 
suppressed dephasing at zero detuning is found below a 
critical value of the Luttinger parameter, which controls 
the correlation decay exponent in the quantum liquid. 
The transition is generated by an orthogonality catastro- 
phe in the coupling to critical one dimensional fluctua- 
tions. We showed that such a non-equilibrium transition 
can be observed in a system of attractive fermions in the 
BCS to BEC cross-over. 
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DIP, ISF, and the US-Israel BSF. 
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